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We investigate the eflective action of the Polyakov loop with spatial variations. We expand the 
effective action not in powers of derivatives or momenta, but in powers of variational amplitudes. At 
one-loop order the results suggest that the instability towards the confining vacuum may be caused 
by the variational terms. 
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(3 . I. INTRODUCTION 

o : 

' It is believed that the vacuum of the SU(iV) Yang-Mills theories should be in the colour-confined phase at zero or 

\ low temperature, while at sufficiently high temperature the colour-deconfincd phase would emerge where perturbative 

^ • calculations do well or at least somewhat make sense. 

O ■ In the SU(A^) pure Yang-Mills theories at finite temperature the centre symmetry is known to be deeply connected 

, with the criteria of colour confinement ||l|,||j3| . The order parameter of the symmetry is given by the expectation value 

QQ of the Polyakov loop operator, 



> ; I 'i'^^^^^)) = # / d'x( ^trTexp ( -i j AxiAi{xi,x) ) ) . (1) 

O 




/3 



Here [3 — 1/T is the inverse temperature, T is the time ordering operator and tr denotes the trace over the colour 
index in the fundamental representation. One can readily relate it to the free energy Fq in the presence of a single 
T-H , quark with infinitely heavy mass: f2 — exp{—/3Fq). Therefore, the confined phase {Fq — oo) corresponds to f2 = 0, 
■ whereas ft ^ in the deconfined phase. If the vacuum is centre symmetric, that is to say fl — Lfl (L G Z{N)), it 
follows that = 0, i.e. that the realization of the centre symmetry is accompanied by the confining vacuum. 
^ The effective potential of the Polyakov loop has been calculated both in the perturbative approach UJ|,§,@,|| and in 

C^H, the strong-coupling expansion on the lattice . The perturbative results have revealed that the centre symmetry 
^ ^ is broken spontaneously due to the appearance of the Debye mass at high temperature. It was also proven that any 
perturbative calculation inevitably leads to f2 = 1, that is the perturbative vacuum pl| ], though the proof is based 
on the suspicious assumption that the Haar measure could be harmlessly negligible. Recently, however, Bohacik 
and Engelhardt and Reinhardt Q have pointed out in the Polyakov gauge (static gauge d^A^ = with an 
additional condition that A4 should be diagonal) that the effective action, which includes not only potential terms 
but also kinetic terms, has a singularity right at the perturbative vacuum. The singularity stems from the Haar 
H [ measure, whose meaning is still a little confused [ p^ , p^Jl6| . This fact seems to suggest that the Haar measure or some 
■ - - 1 configurations with spatial variations such as calorons (multi-instantons at finite temperature), magnetic monopoles 
and domain walls might be relevant to the dynamics of the confinement-deconfinement mechanism. 

These former results, however, are somewhat ambiguous with respect to the two points, i.e. renormalization and 
gauge invariance. Because of the singular character of the Polyakov gauge, it is likely that gauge spurious objects might 
arise from incomplete renormalization. Thus instead of using the singular Polyakov gauge we prefer the covariant 
background field gauge in our calculation as in at the cost of the renormalization necessary for the Polyakov 
loop. 

One might think that the kinetic energy must be dominated by the tree term and that the loop corrections should 
be much smaller than the tree term so as to make sense of the perturbative calculation. However, the longitudinal 
gluon propagator, for instance, acquires the Debye mass from the contribution of the loop corrections and then in 
the soft momentum region the tree term ~ (dA)^ proportional to squared momentum can be certainly negligible 
compared with the finite Debye mass. Of course, when the momentum becomes larger the tree term grows to be 
dominant eventually. Consequently, one has to take account of loop corrections until the momentum exceeds the 



*E-mail: fuku@ntl.c.u-tokyo.ac.jp 
^E-mail: ohta@ntl.c.u-tokyo. ac.jp 



1 



order of T. Because we are interested in the momentum region around T, the effective action should be expanded in 
powers of variational amplitudes instead of derivatives or momenta 13 



The outline of this paper is as follows. We define the effective action of the Polyakov loop in Sec. ||. Using the same 
technique as the derivation of the gauge dependence identity, we show that the effective action of the Polyakov loop 
is generally gauge invariant even if the classical fields do not satisfy the on-shell conditions (the equations of motion) . 



Sec. ^II| is divided into three subsections. Firstly, we expand the one- loop effective action to quadratic variational 
terms. The result suggests possible instability for the perturbative vacuum. Secondly, quartic terms are found to 
stabilize the slipped vacuum. Finally we remark upon the zero temporal modes and the resultant cubic term in the 



potential. In Sec. iV the obtained effective action is examined numerically to give rise to the vacuum instability. We 
conclude our discussion and comment upon problems inherent in our result in Sec. 0. 

II. DEFINITION OF THE EFFECTIVE ACTION 

The effective action in terms of gauge fields is defined in a field-theoretic manner by the equation 

exp {T[A]) = J VAVcDcexp {S[c, c. A] + Ssrc) (2) 

where the action S is the sum of the classical action Sd, the gauge- fixing action Sg[ and the ghost action ^ghostj 
namely, 

S[c, c. A] = Sci[A] + Sgi[A] + S-ghostic, c. A] 

" 4 f 1 1 (5 J-"" j,1 

and the source action S'srcj which only removes the one-particle reducible diagrams, is written as 

Ss.c = f d^xj,4A; - A^^) j^, ^ (4) 

Here ^ is the gauge parameter and J-'"' is to be chosen appropriately to fix a gauge. The upper and lower indices 
of colour should be distinguished, as we employ the ladder basis later (the Cartan metric in the ladder basis is not 
diagonal). For the purpose of defining the effective action of the Polyakov loop, we replace a part of the degrees of 
freedom, that is described by the elementary fields, with the (extended) Polyakov loops defined as 
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f]"(a;) = — trjrexp |^-i J dxiAiixi.x) ) } (5) 

where n runs over l,2,...,iV — 1 corresponding to the Cartan subalgebra (the neutral gluons), or alternatively 
1, 2, . . . , iV with a constraint to force the untraccd Polyakov loops to be Hermitian. Because the Polyakov loops are 
static and written only by the gauge potential A4, it is natural to transfer the variables from A2{x) to fl"{x), where 
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A2{x) = - / dxiA2{xi,x). (6) 
P Jo 



Then the source action (|j) can be replaced by 

f/3 



Sarr / d 



i*x{j^a(A;-A;i)+X„(0"-n")} (7) 



-^f^a^ tit; — I / da;4J4„(a;4,a;) = 



K„ 



ST [A, n] 
sr[A, n] 








where the source J^a is almost the same as that in the source action (^ except that no more A^ arises out of it. 
The expression in the brackets represents this constraint on J^^. The source for A^ is not J^a but Kn, in which the 
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Polyakov loops are regarded as the dynamical variables in place of A". Then setting A = we define the effective 
action of the Polyakov loops that is denoted as T[Cl]. 

In the same manner as the derivation of the gauge dependence identity [0 , it is easy to show that the effective 
action defined above is completely independent of a gauge choice of cither or ^. The infinitesimal change of the 
gauge-fixing condition, 

[A] J'" [A] + bT" {A\ + (8) 

induces the variation only for the gauge-fixing action as the following: 



^ 1 



where Q is the full propagator of the ghost fields defined by the equation of motion, 

'^Sj^^{D^)\{y)g\{v ~z) = -5\5{x - z). (10) 
It is the essential point that (55gt can be compensated by an appropriate gauge transformation, 

SAlix) - {D^)%{x) r d'^yg'^ix - y)5'T^{y). (11) 



The only terms varied under the subsequent transformations (^ and (^jj) are the source terms, to say more precisely, 
fields coupled with the source, provided that the integration measure contains the factor I/VJ [|l^. In the present 
case the source terms are given by (0) among which one term with J^a vanishes because we have set ^ = 0. In the 
other terms with if„, the coupled quantum fields transformed by (11) are the Polyakov loops 17". They are definitely 
invariant for any periodic gauge transformation. Therefore, it is obvious that the effective action of the Polyakov loops 
does not depend on a gauge choice at all due to the gauge invariance of the Polyakov loops even when the classical 
fields do not satisfy the stationary condition. 

For the purpose of the perturbative calculation it is convenient to expand the source term (^ in powers of the 
elementary gauge fields. We parameterize the classical Polyakov loops as 

1 1 ^ 

17" = _tre-™'3° = — V e"'"'^''' (12) 
N N ^ ^ ' 



with the Hermitian constraint. 



N 



N 



J2a, = 0. (13) 



Later on we will make use of the background field method with the background fields a — diag(ai, 02, . . . , ajv) 
introduced in (^). It is important to note that this definition of the background fields contains the renormalization 
of the Polyakov loops from the beginning Q. The gauge fields are decomposed into the quantum parts and the 
background parts respectively, 

A^,{x) = A^,ix) + (5^4a(a;)- (14) 
Using equations (|), (0), (|l|) and (|l4|), we can expand the source term as follows: 

r/3 
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with the notation A^i_j) for the component of the matrix Ai divided by a/2, that is, an off-diagonal 

component of the ladder basis. In the last line is defined by in which A is substituted for A. If we proceed 
to the n-loop calculation of the effective action T, we have to use the {n — l)-loop result for F into ST/Sa and at the 
same time the neglected terms 0{A'^) are also necessary. 

In the final expression of ( [l5| ) the first term simply removes all the one-particle reducible diagrams as in the case 
for the elementary gauge fields written in (^). The second term corresponds to the renormalization of the Polyakov 
loop, which is not the renormalization of the ultraviolet divergence and is necessary for the two-loop or higher-order 
calculations of the effective potential . We emphasize that any effect of the renormalization of the Polyakov loop 
automatically comes out in our formalism because we began with the appropriate definition of the effective action. 
Of course, it is also possible to understand the second term from the point of view of The background field a is 
renormalized to a + Sa and then the effective action can be expanded as r[a-|-5a] = T[a] +ST/6a-Sa, the second term 
of which is the same form as the second term in (|l^) . In the present case of the effective action the situation is slightly 
involved. The zero-loop order of the kinetic term (tree term) may survive so that the effective action of the one-loop 
order can also be affected by the renormalization of the Polyakov loop with spatial variations. However, although 
finite corrections may be certainly generated, they can be absorbed into the arbitrariness of the renormalization of 
ultraviolet divergence and thus we do not have to bother with the renormalization of the Polyakov loop up to one- loop 
order at all. In the next section we explore this point further. 



III. PERTURBATIVE EVALUATION OF THE EFFECTIVE ACTION 



The background field method has an advantage in that the effective action preserves the manifest gauge invariance 
[ p8[ . The gauge fields are decomposed into 

^m=-^m+^Bm (16) 

where and denote the quantum and the background part, respectively. Introducing the covariant derivative 
in terms of the background fields by Db^ = 9^ + i[^B^t, ], one can write the gauge- fixing condition as 

Db^a;-lu^^o. (17) 

Multiplying a Gaussian weight and integrating out, we obtain the Lagrangian, 

-Ah^DB^A^aAlAl + f^thecA^AlAlAl] + [{dI)\ + {DB^^VbdK] ■ 

From the quadratic part of this Lagrangian, we can readily integrate the quantum fields and in the Feynman 
gauge ^ = 1 to give the one-loop effective action, 
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(19) 



where tr is the trace of the spatial and colour indices, while tr' includes the trace over the Lorentz indices in addition. 
We set the background field as (|lj) and a is decomposed into 

a{x) ^ ^{C + e{x)) (20) 

where C and e{x) represent the constant and the variational part of the background field, respectively. They are 
allowed to be traceless independently, i.e. 

N N 

trC = ^ C, = tre(a;) = ^ e.,{x) = 0. (21) 

i=l 1=1 

We divided them by (3 in order to make C and e{x) dimensionless. In the ladder basis j^,^ the effective action ( pj| ) 
is attained most simply as the following: 

r[C, e] = -g^j^ Yl ~ ^tr' In 01^)5^,^ + 2i {(9^(5^4 - (5p49,.)ey } 

' ^ i>3 

+2\eij5^udi(^i,j) - [eijf'S^iy] + trln + 2iey94(ij) - (cy )^| (22) 

where we rescaled all the coordinate variables by (3 so as to let them be dimensionless. From now on, we consider all 
the coordinates, momenta and fields as dimensionless multiplied appropriately by (3. Here we denote the dimensionless 
derivatives with the background fields as 

Cij — C'i ^ Cj ^ij — £i £j- (24) 



Our purpose is to expand the effective action (^2|) in powers of e. In terms of Feynman's diagrammatical method, the 
terms with e and in the action (^2|) are interpreted as the three- and four-point interaction vertices. The quadratic 
derivative terms correspond to the free propagator. 

The one-loop effective potential shown diagrammatically in Fig. |^ has already been argued in |^ with the result, 

i>j 



mod2 



where n in the summation is the Matsubara frequency, which appears only through the momentum variables P(ij) — 
{27rn + Cij,p). 

A. Up to quadratic order 

The first correction up to order is represented diagrammatically in Fig. |^. Introducing the Fourier transform of 
the variational fields, 

= / ^^^^(l)^^'"'" (26) 



we can easily write the tree term as 



r.ee[e]=-^E/A^|..(.)|V. (27) 



5 



The procedure to calculate the loop corrections is almost the same as that in [Q except for the absence of the 
condition that the momentum is sufficiently small. This condition is really the case for the Z(7V) domain wall but in 
contrast with it our results should contain all orders of the derivative expansion so that it would not take the form of 
the ordinary kinetic term proportional to momentum squared. The renormalized contribution to the effective action 
is calculated as 



i.j ^ ' n 



(27r)3 I + g)2 r pI 



V ^ f , , ,,2 f 11 2, 1 

/ (2', 

i>'j 



where we used the 4- vector notation q = (0, q). M is a renormalization point. The first term in brackets in ( [2q ) 
comes from the vacuum part and the second one comes from the matter part which vanishes as T ^ 0. 
The function f{q, C) is defined by 

2 r._L,^Y, (29) 



f{q,C) = -l dp^p + g( l-^)ln 
1 1 



2p + q 



^(^' ^-■) = ^Fhc— T + (30) 

whose behaviour is shown in Fig. ^ N{p, C) is the distribution function which is reduced to the usual Bose distribution 
function at C = and to the negative Fermi distribution function at C = tt. This exchange of the statistical character 
agrees with the centre-symmetric treatment proposed recently in which the charged gluon must obey the antiperiodic 
boundary condition jl^] (of course, the origin is actually identical). If the coefficient of N(j),C) is positive definite 
and also sufficiently large, the minimum of the effective action will shift into C ~ tt, i.e. the confining vacuum. In our 
formalism the result of |l^,^ can be interpreted as that non-zero value of |eij(q')| causes Cij ^ tt as the temperature 
is lowered. However, as we will discuss in the subsequent sections, our consequence suggests another scenario: dj 
remains zero and only |eij(q)| makes the Polyakov loop decrease. Furthermore, the behaviour in the vicinity of C = tt 
will prove to be so infrared singular for spatial variations that it is out of reach of the perturbative regime. 

We have to explain now that the renormalization of the Polyakov loop is actually harmless to the above result, as 
mentioned before. The second term of (^5|) corresponds to the first correction of the Taylor expansion in the variation 
from the renormalization of the Polyakov loop, that is, 5e. It is possible to calculate Se from ([l5| ) at the lowest order 
of e to convince ourselves that Se is independent of Cy so that the correction to the effective action takes the form of 
(constant) xg^ I ey(q)p for the quadratic order of e. This can be absorbed by the redefinition of M. Consequently, we 
can safely ignore any effect of the renormalization of the Polyakov loop up to one-loop order. 

The above result can be compacted into 

^ -^Y.J i^i'.Mi' + ln^.c.,}} (31) 

using the running coupling constant, 

^ ~ l + {llNg^/A8TT^)HqyAP) ^ llNlniqyp^Al) ~ nN\n{qt) ^ ^ 

where t denotes the dimensionless temperature in units of Aq, the QCD scale parameter. 
Asymptotic values of f{q, C) are easily obtained as 

/(g^0,C) = ^~4C-H— (33) 
/(q->oo,C) = --/(g-^0,C) = -27r + 6C (34) 

Z TT 

It is obvious that ( ^3|) at C = provides the correct perturbative result of the Debye mass which stabilizes the 
perturbative vacuum. The behaviour of (|4|) that /(g, C) turns to be negative as q becomes larger is crucial in our 
argument. The perturbative vacuum (Cy = ^ijio) — 0) may become unstable supposing that the decreasing rate of 
f{q,C) overcomes the increasing rate of the tree term (quadratic in momentum). Actually, it is the case when the 
coupling constant is beyond some critical value. 
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Both Cij and are determined by the condition that — r[C, e] (the free energy in thermodynamics) should be 
minimized. From (|3l|) it follows that the vacuum stability for a given e{q) around the perturbative vacuum depends 
on the sign of 

with Q!((7,t) = (7^(g, t)/47r. Of course tij is not independent of each other under the constraint (|l3|). At C = 0, 
however, the coefficient function t) contains no colour dependence (zj dependence) so that a negative h{q, t) is 
sufficient for the occurrence of the vacuum instability. The behaviour of h{q,t) in the case of SU(3) {N = 3) as the 
temperature t varies is shown in Fig. ^. 

Perturbative calculations are not reliable at all for a{q,t) > 1. Looking at Fig. ^ we realize that the momentum 
region with a > 1 indicated by the broken curves grows when we lower the temperature and at the same time the 
possible unstable modes in which h{q, t) < come to spread. Therefore, whether our perturbative evaluation is 
reliable or not near the unstable modes depends on whether the spreading rate of the unstable modes overcomes the 
growing rate of the unreliable momentum region. In fact. Fig. ^ indicates that near the critical temperature where 
unstable modes begin appearing, the relevant modes are in the momentum region for a < 1. To be more precise, at 
the critical temperature — 2.77Ao, the unstable mode is qc = 1.33 and the corresponding critical coupling constant 
is ttc = 1.31/A^. In the case of SU(3), for example, ac = 0.44 which is surely smaller than one. 

However, it is subtle to decide from the critical value of the coupling constant whether or not the perturbative 
calculation works well at that point. Nevertheless it is most likely that this instability actually occurs because we 
can comprehend it in an intuitive manner as follows. Under the effect of the temperature, real (on-shell) particles 
are excited thermally. The average distance between the thermal particles is of order of 1/T and their distribution is 
incoherent. Owing to this incoherency the vacuum tends to be stable against spatial variations with small momen- 
tum or equivalently with larger-scale variations compared with the average distance between thermal particles. For 
variations with momentum larger than ~ T, on the other hand, the vacuum tends toward instability as the overlap of 
the vacuum configuration with the distribution of thermal particles increases. If the strength of the interaction, that 
is to say, the tendency to overlap the configuration is large enough, the vacuum may become unstable. 

Once unstable modes emerge, the variational fields e^j go on increasing to minimize the free energy unless the higher 
terms are taken into account. In order to investigate the property of the instability more closely we need higher-order 
terms of e^j. This is the subject dealt with in the next subsection. 



B. Up to quartic order 

Corresponding diagrams to be evaluated are displayed in Fig. ^j. The result obtained after tedious calculations has 
singularities at C = 0, whereas the quartic contributions are ultraviolet and infrared finite as long as C 7^ 0. The 
reason for the presence of singularities is quite obvious: since we expanded the effective action in powers of e/(C or T), 
our expansion is not valid at C = for the temporal zero- mode for which T is absent. The centrifugal barrier at C = 
observed in p^ , p^ seems to arise from these spurious singularities. Therefore, we must manage infrared divergence 
by resumming zero modes, as in the case of ring resummation for infrared divergence. We remark upon this problem 
in the next subsection. 

In this subsection we evaluate the variational term of quartic order from which all the temporal zero-mode terms 
are subtracted to be resummed afterwards. Once we get rid of the zero mode it can be expected that the result will 
be almost independent of q in the interested region q ^ 1.3, because the expansion in terms of q becomes a power 
series of 

'^"^ -~ (36) 



2 



(27rn)2+_p2 {2Tm)^+p 

where p is the momentum to be integrated and n is the Matsubara frequency. Although in the case of n = they 
are singular near p '--^ 0, we can regard them as small enough as far as n 7^ and q ^ 1.3. Of course, it is possible to 
compute the quartic correction directly from the corresponding diagrams. We confirmed numerically that the result, 
which is too intricate an expression to write down here, hardly depends on q around the unstable region. 

Now that the result proved to be almost g-independent we can evaluate the quartic correction at the limit of g — s- 0. 

^ [^'^] -^-psl^J (2^(2;^(2^'^^'^'^i)'''^('^2)e^,(g3 - gi)ey(-92 - Qs) 
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(2^ 



E 



4p4 



P7 ■^ 



Vy 



(27r3) 

-oo+ir; (27'') 



V_ 



E 



1>J 



(27r)3 



d3gi d3g2 d3g3 
(27r)3 (27r)3 (27r)3 









d3p 




(27r)3 




8pl 1 ■ 






( 


d3p 1 


93 


(27r)3 12 dp^ 



2p4 



1 



1 



1 



1 



(-,-ip4-iCi. 
1 



1 



-iCi. 



127r2 



(37) 



The first, second and third term in the square bracket corresponds to Fig. || (a) — (c), respectively. The result (|37| ) is 
nothing but the coefficient of in the Weiss potential (^5|) . Since the Weiss potential is a quartic polynomial, we can 
realize without any calculation that terms of order higher than quartic are never relevant under the soft momentum 
limit if subtracted by their zero mode. 



C. Zero- mode resummation 



One may wonder why we did not take account of the cubic terms of e in spite of the fact that the Weiss potential 
(^) has the cubic term . If we compute them perturbatively as in the same manner for quadratic or quartic order, 
however, they are always zero: the integrand for the cubic order is proportional to Pi{ij) so that the summation 
over the Matsubara frequency and the colour indices leads to zero. Then where does the cubic term in the potential 
( p5| ) come from? The answer is the zero-mode resummation, which is quite similar to the e3-correction of the QED 
partition function ||l9t| . 

Here we will perform the zero-mode resummation only for the most dangerous terms that are divergent in the 
infrared limit. The total momentum insertion from the variational fields e is necessarily zero for one blob (self-energy 
part) so that each propagator from one blob to the other possesses the same momentum to bring about the strongest 
infrared singularities, as is represented in Fig. ^. 

Taking away all the terms with n =/= and subtracting quadratic terms (infrared finite terms), we reduce the effective 
action (E2[) into 



tr' In 

4 



Sni, 



2e. 



'^5 - 



a2 f)2 "^M 



Sui, + 



4{(5p(5A4 - Sf_,4dx)eij} {{dxS„4 - S\4di^)eij} 



f)2 



id) 



B^ 



-— tr In 
2 



2f2. 



+ 



(e?- -terms) 



f)2 



.2 



"{^.3) 



(38) 



where any term of is absent because all the diagrams possessing such vertices are to be resummed individually to 
contribute either to the order of or to less singular terms. In order to evaluate ( |3^ ) we perform calculations in the 
momentum space. e{x) is sum of all the modes e(q), among which we leave only the most singular combination, as 
shown in Fig. ||. 

To take the trace for the Lorentz indices, it is useful to decompose them into three projective components, namely 



4/33 

1 



(27r)3 2^ 



In 



Sup + 



_f_q 

{2tt 




'/jiy = Si^Siij. Then ( |38D can be rewritten as 
2 



46*2. 



VI3) PhM^,3)+^?) 



i2n) 



q 



3l^t3 
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1 



V 



W) 




— (e.^j-terms) 



P%,3) P%jM^,j)+<lf 



(27r) 



2 4(C^-+a 



(39) 



We restrict ourselves to considering the case of Cy = because no infrared singularity emerges otherwise. For an 
isolated mode q, the zero-mode action (|39|) is simplified as 



'2/33 (27r)3 



ln<^ 1 



4q2 



(40) 



It is easy to reproduce the cubic term in the Weiss potential ( |25| ) with q = 0. In this limiting case ( ^ ) is reduced to 



V 



y- 



dp- 



/?3 



37r 



(41) 



This is certainly identical to the cubic term in the Weiss potential. 

If the relevant mode q is non-zero then the situation changes drastically. For sufficiently small values of je^ | ^ q 
the main contribution of the integration over -p stems from singularities near p ^ ^ I- Then a similar calculation 
as above leads to 



'^53^2^ JO 



dp \ In 1 



V ^ 1 



dp- 



(42) 



The integration in the last line is defined in the sense of Cauchy's principal value. The result (^2|) is, in fact, what 
is expected: the e3_term is generated by the most severe infrared singularities. Once q takes a non-vanishing value, 
infrared singularities within the blobs (not the propagators, see Fig. ^) are loosened. 

The result obtained up to now is valid only for le^jl ^ q. In the next section, we want to deal with not so small e 
around q^, beyond which the perturbative evaluation will lose its validity as shown in Fig. 0. If the contribution from 
the zero-mode resummation comes to be opposite to that from the quartic term for some marginal value of e, the 
zero-mode resummation may cause qualitative changes: the absolute minimum of our effective action may jump with 
a finite gap. In fact, however, the zero-mode resummation rapidly grows as q increases from zero so that it becomes 
positive contribution (the same sign as the quartic term) in the momentum region of interest. No qualitative picture 
can be affected by the zero-mode resummation. Thus for the simplicity we will entirely ignore this contribution 
in the next section. Although this neglection seems too rough, closer investigation of the zero-mode resummation 
reveals that it is not so crude. Numerical cancellation between the resummed contribution performed here and the 
next-to-leading contribution occurs. 

The expansion singularities of e/C resummed in this subsection have much to do with the famous problem of Lindc's 
infrared catastrophe [ po[ . Since Fig. ^ is the same as the diagram for the gluon vacuum polarization, the two- or 
higher-loop calculation should break down and non-perturbative effects are no longer under control perturbatively. 
The infrared singularities would be regulated and perturbative expansion would work well only provided that the 
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constant background field C is of order of (7° or at least as is discussed in g. It is the reason why the hard thermal 
loops (HTLs) are not necessary for the evaluation of the effective potential. The naive order counting by the coupling 
constant will work due to the presence of the background. Since the stationary value of the background is zero, in 
fact, there seems to be a logical contradiction in this argument. In our formalism the variational background field e 
instead of C takes a non-zero value, as seen below. Linde's problem is regulated here into the expansion singularities 
of e/C that can be resummed in principle. Therefore, we can consider that our perturbative expansion of the order of 
g should make sense (this does not say nothing about the reliability of the perturbative evaluation) and expect that 
a non-perturbative contribution will be suppressed once instability occurs. 



IV. NUMERICAL ANALYSIS 



The free energy density obtained above is 

/33 



c/[c, e] = - ^ (rw [C] + r(2) [c, 6] + [6]) 



12 \ TT 



1 f d'\i d^q2 d^qs / n . n . n / n 

12^ J (2.)3 (2^)3 (2^)3 ^-.(gi)^-.(9^)^-.(g3 - gi)^.(-g2 - 93) 



(43) 



We can acquire the equations of motion from the stationary conditions for the free energy, i.e. SU/6Cij = SU / 5eij{q) = 
and, in principle, solve them self-consistently. It is, however, quite involved because of the mode couplings induced 
by the nonlinear term e**. To simplify the analysis and investigate the situation in a comprehensible manner we assume 
that e^(a;) is such a function that the variance from the spatial mean value is small enough to be negligible, that is, 

^ J d'x |4(a.) - 1 1 d3x4(a.)|' ^ 0. (44) 
Then the nonlinear term takes the form, 

This assumption is correct provided that e{x) is a step-like function or has no sharp peak. 

For simplicity we shall initially discuss the SU(2) case where £2 = ~t\ = e/2 and C2 — ~C\ — C /2. From the 
equations of motion one can easily obtain 

'^''^ 'e(q)|2 = -6^2^(C,t), (46) 



(27r) 

where h{C, t) is defined as 

with (7o(C, t) determined in such a way that q'^ /Na{q, t) + f{q, C) should take its minimum value to minimize U[C, e]. 
Then the free energy density is written only in terms of C as 

U[C] = U^''\C]~i^''{h{C,t)Y, (48) 

where JJ^"^ [C] is the Weiss potential (density) , 

^'"'[^i^f?)^ {(7) -F- <^») 

12 Jraod2 I V / mod2 J 

The order parameter (|l|) takes the following form under the approximation (E3) and using the stationary value (E^) , 
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n = exp(^^h{C,t)^cos^. (50) 

Minimizing (^sj) leads to C = so that ft varies only through the change of e determined by (^ ) at C = 0. The 
numerical behaviour of the order parameter as a function of t is displayed in Fig. |7[ 

It is straightforward to proceed to the case of SU(3) where the independent variables are ei, £2, Ci and C2. The 
traceless constraint (^3|) gives €3 — — ei — £2 and C3 = — Ci — C2. Numerical investigation tells us that the minimum 
lies either at £1 — £2 or at ci = — 2£2 (£2 — — 2£i). Whichever of them we may choose, they are actually equivalent. 
We take the former case, namely, £1 = £2 = £• The stationary value in this case corresponding to (Eq) is 



2 

HQ)f - -y {hC^Ci + C2,<) + h{Ci + 2C2, t)} , (51) 



(27r)3 

and the free energy density corresponding to ( ^ ) takes the form of 

U[C] = C/(°) [C2 - Ci] + C/(") [2Ci + C2] + C/(°) [Ci + 2C2] 

■^7^2 r _ , -.2 

-— {/i(2Ci + C2,i) + /i(Ci + 2C2,t)} (52) 

whose minimum turns out to be at Ci = C2 = 0. The behaviour of the order parameter shown in Fig. ^is almost the 
same as in the SU(2) case. 

It should be noted that we have discovered the possibility of instability but we cannot reach the phase transition. 
The phase transition is defined by the critical point above which spontaneous magnetization (order) appears. Any 
critical characteristics, such as the transition temperature, critical exponents and so on, are all defined at the very 
critical point. Thus all we can say is that Fig. ^ suggests some continuous instability. If we want to investigate the 
critical point, we must start from the confining vacuum as the strong-coupling limit on the lattice, for example, or 
some non-perturbative formalism [p^ . 

As mentioned above, our result of the free energy density ( p3| ) cannot describe any discontinuous transition which 
is believed to be the case for SU(3). The reason for this is obvious: equation ( ^3|) does not possess any terms of e^. 
Since the centre group Z(3) allows terms of fl^ only in the case of SU(3), terms are expected to be produced from 
some origins peculiar to SU(3). 



Another possibility to lead to the discontinuity is the instability near C ^ n. Even if (43) has no odd terms of 
£ it is a complicated nonlinear function in terms of C so that the discontinuous transition in terms of C cannot be 
ruled out. Actually, Fig. ^ shows that /(g, C = n) becomes negative in the soft momentum region, which may cause 
h{q,t) to be negative. Although this instability seems to be similar to the case of C ~ 0, the relevant modes are 
entirely different. On one hand, in the case of C ~ 0, the unstable modes are around g ~ 1.3 that is comparable to 
the temperature. In the case of C ~ tt, on the other hand, the instability is enhanced as the momentum goes to the 
soft limit q ^ 0, where the running coupling constant at the one-loop order does not make sense. Our result that ( |5^ ) 
has the absolute minimum at Ci = C2 = would be valid provided that the minimum value of h(C — 7r,i) were at 
most f{q — > 0, C = 7r)/47r = —1/6 (we have used this assumption implicitly in searching for the minimum of ( ^ ) and 
(^)). However, the behaviour of h{C,t) obtained here is only up to one-loop order and thus we have no idea what 
really happens near C ~ tt. All we can state is that C ^ tt is certainly unstable and that the relevant momentum is so 
infrared that any perturbative calculations would not do well. It is therefore probable that the position of the absolute 
minimum jumps from C = to tt in the SU(3) case if the free energy density is attained in a non-perturbative way, 
which is beyond the scope of this paper. 

The situation is clear in Fig. |^ which shows the behaviour of (^8|). The instability discovered here is for the modes 
near C = at lower temperature. Around C ~ tt the relevant modes are so infrared that non-perturbative effects 
turn out to be of importance, even if the temperature is sufficiently high. It might seem curious that in Fig. ^ the 
potential value at i = 3.0 is lower than that at t = 1.6 near C = tt; this is because we have searched for the most 
unstable mode qo within a range larger than l/t in order to avoid spurious negative values of the running coupling 
constant. When the temperature becomes larger, the lower bound of the momentum range to be searched decreases 
and thus the potential value also decreases. 

V. SUMMARY AND CONCLUSION 

We investigate the effective action of the Polyakov loop defined in a field-theoretic manner. The gauge invariance 
of the effective action can be proven by the same strategy as deriving the gauge dependence identities. Owing to the 
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appropriate definition of the effective action any effect of the rcnormahzation of the Polyakov loop that is necessary 
for higher-order calculations is automatically generated through non-local interaction vertices. In order to evaluate 
the effective action we choose the background field gauge where the background gauge symmetry remains manifestly. 
In this gauge the most essential properties to attain reliable results, that is, the gauge invariance and the unambiguous 
renormalization, are satisfied simultaneously. 

In the actual calculation we expand the effective action in powers of the variational amplitudes e(x) instead of the 
derivatives. The quadratic term of the variational field shows the possibility of the vacuum instability for marginal 
modes around C = or for the constant background C ~ tt. The reason why the vacuum may become unstable 
is comprehensible from our resultant expression: as for the former case (C = 0), the variational backgroimd is 
usually absent because such a configuration costs more kinetic energy. In the thermal circumstance, however, on-shell 
particles are excited to obey the thermal distribution function, in addition to the quantum excitation. If the mode of 
the variational background becomes comparable with the inter-distance of the thermal particles, i.e. ^ 1/T and if the 
interaction is sufficiently strong, then the vacuum can decay into non-vanishing variational background. As for the 
latter case (C ^ tt) the statistical character of the transverse gluon changes into that of the Fermi particle, that is to 
say, the transverse gluon obeys the antiperiodic boundary condition. Moreover the distribution number of particles 
turns to be negative. The serious problem we encounter is the magnitude of the coupling constant. It depends on 
the smallness of coupling constant whether our effective action at one- loop order is reliable or not. In the SU(3) case 
the critical coupling constant is etc = 0.44. Particularly at finite temperature this value is too large to regard it as 
small. It is well known that the higher-order corrections to the free energy up to order (this is the upper limit 
of Linde's problem) is so oscillating that the perturbative series does not seem to converge at all. Nevertheless we 
believe that the physical picture may be refiected qualitatively in the one-loop effective action. It is an interesting 
and essential point that the kinetic term of the effective action even at one- loop order can bring about instability at 
the perturbative vacuum, which is out of reach of the conventional effective potential. The effective action at one-loop 
order already has a non-trivial physical content. 

Once we pursue the policy of studying the physics contained in the one-loop effective action or believe in the 
perturbative evaluation despite the large coupling value, we can proceed to expand the effective action to quartic 
order. A brief discussion reveals that the quartic term is independent of modes provided that the temporal zero mode 
is subtracted. Looking over the Weiss potential we can see that any terms higher than quartic order come from the 
temporal zero mode. The zero-mode resummation is performed, where we make clear the origin of the cubic term in 
the Weiss potential. 

The effective action density determines the stationary vacuum for the variational background e{x) and the constant 
background C. The minimum lies at C = so that the Polyakov loop decreases only by no n- vanishing e provided 
that the instability near C = tt is harmless. Since we know that in the SU(2) case the transition is continuous (second 
order), the instability found in this paper is likely to be relevant to the real physics. In the SU(3) case, on the other 
hand, the discontinuous transition implies the possibility that the infrared instability around C = tt comes to play an 
essential role. 

It is necessary, of course, that the qualitative argument discussed here is justified by some reliable method other 
than perturbative evaluations. Some methods, such as the thermal renormalization group and infrared effective theory 
successfully predict the non-pcrturbative contribution for the infrared region. The application of these methods to 
estimate the effective action is quite difficult due to the presence of the non-local intricate vertices corresponding to 
the renormalization of the Polyakov loop. That is the future work to resolve. 
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FIG. 1. Diagrams for the one-loop effective potential. The curly curves and the broken curves represent the gauge fields and 
the ghost fields, respectively. 

FIG. 2. Diagrams up to quadratic order. The doubled curly curves denote the variational fields. 
FIG. 3. The behaviour of the function f{q, C) with various values of C . 
FIG. 4. The behaviour of the function h(q, t) for various temperature t. Broken curves represent the momentum region 
where a > 1. 

FIG. 5. The variational corrections of quartic order. 
FIG. 6. The ring diagram : each pair of doubled curly curves has the opposite momentum so that every single curly curve 
possesses the same momentum (inserted momentum is always zero). 

FIG. 7. The behaviour of the order parameter as a function of the temperature t for SU(2) and SU(3). 
FIG. 8. The free energy density only in terms of C: The broken curves are in the unreliable region a > 1. 
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